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Abstract. Carrier-induced ferromagnetism has been observed in several (III,Mn)V 
semiconductors. We review the theoretical picture of these ferromagnetic semiconduc- 
tors that emerges from a model with kinetic-exchange coupling between localized Mn 
spins and valence-band carriers. We discuss the applicability of this model, the validity 
of a mean-field approximation for its interaction term widely used in the literature, and 
validity limits for the simpler RKKY model in which only Mn spins appear explicitly. 
Our conclusions are based in part on our analysis of the dependence of the system's 
elementary spin excitations on carrier density and exchange-coupling strength. The 
analogy between this system and spin-model ferrimagnets is explored. Finally, we list 
several extensions of this model that can be important in realistic modeling of specific 
materials. 



1 Introduction 

Semiconductors and ferromagnets play complementary roles in current infor- 
mation technology. On the one hand, low carrier densities in semiconductors 
facilitate modulation of electronic transport properties by doping or external 
gates. Information processing is, therefore, most commonly based on semicon- 
ductor devices. On the other hand, because of long-range order in ferromagnets 
only small magnetic fields are necessary to reorient large magnetic moments and 
induce large changes in magnetic and transport properties. This explains the 
use of ferromagnetic metals in storage devices. The application of giant mag- 
netoresistance (GMR) and tunneling magnetoresistance (TMR) effects to read 
magnetically stored information is the poster child of spin electronics in ferro- 
magnetic metals. The prospect of synergism 0J^] between electronic and mag- 
netic manipulation of transport properties has raised interest in spin electronics 
in semiconductors. In a (III,Mn)V ferromagnetic semiconductor charge and spin 
degrees of freedom could, it is hoped, be manipulated on the same footing in 
basically the same material used now in highly successful electronic devices. 

Two preconditions for spin electronics in semiconductors are, first, the gen- 
eration of spin-polarized carriers and, second, their spin-coherent transport in 
the material. While spin-coherent transport over several micrometers has been 
demonstrated HU^], spin injection from ferromagnetic metals into semiconduc- 
tors has proven to be problematic || . More promising is the use of ferromagnetic 
semiconductors, in which not only charges but also magnetic moments are in- 
troduced by doping pip). 
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1.1 Diluted magnetic II- VI and III-V semiconductors 

Semiconductors have rich optical and magnetic properties when Mn or other 
magnetic ions are doped into the material For a long time, the investi- 

gation of these so-called diluted magnetic semiconductors (DMS) concentrated 
on II- VI compounds such as CdTe or ZnSe, where the valence of the cation is 
identical to that of Mn. The magnetic interaction among the Mn spins is in 
this case, however, dominated by antiferromagnetic direct exchange, and the 
systems show paramagnetic, antiferromagnetic or spin-glass behavior. To obtain 
ferromagnetism, free carriers must be introduced as we explain below. This can 
be achieved by additional doping fl~j)| , which is, however, rather difficult experi- 
mentally. Only low transition temperatures have been realized by this route so 
far. 

An alternative approach is to dope Mn in III-V semiconductors such as 
GaAs |ll|P,P,p^|P,[i^,P,p^| j2(i|,P H . Ferromagnetic samples with transition 



temperatures T c up to 110 K | fl4|l5|| have been realized using low-temperature 
molecular beam epitaxy (MBE) growth, and it is hoped that still higher transi- 
tion temperatures will be achieved in the future. Due to the difference between 
the valence of Mn 2+ and the cation (e.g. Ga 3+ ), free carriers (holes) are then gen- 
erated in the valence band. The main focus of this article is on the mechanism by 
which these carriers mediate interactions that lead to long-range ferromagnetic 
order among the Mn spins. 



1.2 Model for Mn doped III-V semiconductors 

We use a partially phenomenological approach, starting from a relatively simple 
but soundly based model with effective parameters that can be determined by 
comparison with experiment. We explore the physical properties implied by this 
model. The reliability of the model will ultimately be determined by compari- 
son with future experimental findings. We will start by discussing the simplest 
version of this model; a variety of extensions can be accommodated without any 
fundamental changes. Some of the anticipated extensions will be discussed at 
the end of this article. 

It is believed that Mn 2+ impurities substitute randomly for the Ga 3+ cation 
in the zincblende structure. We assume that the Mn ion, with its half-filled d 
shell, acts as a S — 5/2 local moment. If interactions did not suppress fluctua- 
tions in the Mn ion charge sufficiently, its d-electrons would have to be treated 
as itinerant, and this model would fail. We denote the Mn ion positions by R[. 
Hybridization between the localized Mn d-orbitals and the valence-band orbitals 
can then be represented by an antiferromagnetic exchange interaction. It turns 
out that the net valence-band carrier concentration p is much smaller than the 
Mn ion density Nun ]l5] . This might be due to antisite defects (As ions sitting 
on a Ga lattice site, providing two free electrons per defect) which are com- 
mon in low-temperature MBE grown materials. As will become clear, we believe 
that this property of the as-grown materials is actually essential to the robust- 
ness of their ferromagnetism. These considerations suggest the following model 
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Hamiltonian: 

H = H a + J pd Jd 3 r S(r) ■ s(r), (1) 

where S(r) = SiS(r — Rj) is the impurity-spin density. In the following, 
we approximate the sum of delta peaks by a smooth function, thus neglecting 
disorder induced by randomness of Mn sites. This approximation is partly jus- 
tified by the carrier to impurity density ratio p/Nmu, since the typical distance 
between Mn ions is smaller than the typical distance between free carriers. The 
itinerant-carrier spin density is expressed in terms of carrier field operators by 
s(r) = | J2crcr< &a(. r ) T <?<?' , P<?'( r ) where r is the vector of Pauli spin matrices. The 
antiferromagnetic exchange interaction between valence-band holes and local 
moments is represented by the exchange integral J p( j > 0, an effective parameter 
which is dependent on the details of the system's atomic length scale physics. 
The contribution Hq includes the valence-band envelope-function Hamiltonian 
[ p3| and, if an external magnetic field B is present, the Zeeman energy, 

H = Jd 3 r^^(r) ^-^-^^(r)-HBB.[ 9 S(r)+g*s(r)]y (2) 

The effective mass, chemical potential, and ^-factor of the itinerant carriers are 
labeled by to*, (j,, and <?*, respectively. To simplify the present discussion we use 
a generic single-band model with quadratic dispersion. 




Fig. 1. Model for Mn doped III-V semiconductors: local magnetic moments (Mn 2+ ) 
with spin S — 5/2 are antiferromagnetically coupled to itinerant carriers (holes) with 
spin s — 1/2. 



We note that the model we use here is identical to those for dense Kondo 
systems, which simplify when the itinerant-carrier density p is much smaller 
than the magnetic ion density Nmu |^|- The fact that p/Nyin *C 1 (while the 
Kondo effect shows up in the opposite regime) in ferromagnetic-semiconductor 
materials is essential to their ferromagnetism. Similar models have been used 
for s — f materials p6| with ferromagnetic exchange and for ferromagnetism 
induced by magnetic ions in nearly ferromagnetic metals such as palladium |27J . 
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Additional realism can be included in the model by using a realistic band model, 
and by accounting for carrier-carrier and carrier-dopant Coulomb interactions. 

1.3 Outline 

The outline of this article is as follows. In Section [2] we point out that neither the 
RKKY picture nor the mean-field approximation to the above kinetic-exchange 
model yield completely satisfactory theories of ferromagnetism in doped DMSs. 
To develop a more complete picture of these ferromagnets we present in Section 
H a recently developed theory |2£||29|,[30| , in which we identify the elementary spin 
excitations (Section |j) of the kinetic- exchange model. We compare the excitation 
spectra that result with those predicted by the RKKY picture (Section |1), by the 
mean-field approach (Section)^), and by a ferrimagnetic lattice model (Section^). 
Section |§| contains an alternative derivation of the spin-wave dispersion, giving 
also some insight in the nature of the excited states. Implications of these results 
for the excitation spectrum for other properties of these materials and extensions 
of the simple model presented above are discussed in Section 0. 

2 Origin of ferromagnetism 
2.1 RKKY interaction 

Simple perturbation theory can be used to describe the familiar itinerant-carrier- 
mediated Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction between local 
magnetic moments. At first glance, this interaction might explain ferromag- 
netism in DMSs, too: each local impurity polarizes the itinerant carriers nearby, 
an adjacent Mn ion experiences this polarization as an effective magnetic field 
and aligns, for small distance between the Mn ions, parallel to the first local 
moment. 

The RKKY theory is, however, a (second-order) perturbation expansion in 
the exchange coupling J p( j, i.e., the picture is applicable as long as the per- 
turbation induced by the Mn spins on the itinerant carriers is small. As we will 
derive below, the proper condition (at zero external magnetic field B) is A <C tp 
where A = Nyi^JpdS is the spin-splitting gap of the itinerant carriers due to 
an average effective magnetic field induced by the Mn ions, and is the Fermi 
energy. This condition is, however, never satisfied in (III,Mn)V ferromagnets, 
partially because (as mentioned in the introduction) the valence-band carrier 
concentration p is usually much smaller than the Mn impurity density A^Mn- 
Instead, the "perturbation" is strong. Valence-band spin splitting comparable 



in size to the Fermi energy has been confirmed by direct experiment |31 . The 
RKKY description does not provide a good starting point to describe the ordered 
state in ferromagnetic DMSs. 

A related drawback of the RKKY picture is that it assumes an instantaneous 
static interaction between the magnetic ions, i.e., the dynamics of the free carriers 
are neglected. As we will see below this dynamics is important to obtain all types 
of elementary spin excitations. 
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Within mean-field theory, however, the spin-splitting of the valence bands will 
vanish as the critical temperature is approached. It turns out that as a result the 
RKKY picture and the mean-field-theory approximation to the kinetic-exchange 
model make identical predictions pl| for the critical temperature. We comment 
later on the accuracy of these critical temperature estimates. 

2.2 Mean-field picture 

The tendency toward ferromagnetism and trends in the observed T c 's have been 
explained using a mean- field approximation p^ , ^3|j34| , |35| , p6| , ^7| , analogous to the 
simple Weiss mean-field theory for lattice spin models. Each local magnetic mo- 
ment is treated as a free spin in an effective external field generated by the mean 
polarization of the free carriers. Similarly, the itinerant carriers see an effective 
field proportional to the local moment density and polarization. At zero temper- 
ature, the total energy is minimized by a state in which all impurity spins are 
oriented in the same direction, which we define as the positive z-axis, and the 
itinerant-carrier spins are aligned in the opposite direction. The spin-splitting 
gap (at zero external magnetic field B) for the free carriers is A = NyinJpdS, 
and the energy gap for an impurity-spin excitation is £} MF = p£J p( j/2, where 
£ = (pi —p-{)/p is the fractional free-carrier spin polarization. 

This picture, however, neglects correlation between local-moment spin config- 
urations and the free-carrier state and, therefore, fails to describe the existence 
of low-energy long-wavelength spin excitations. As we discuss below, correla- 
tion always lowers the energy of collective spin excitations (spin waves) in 
comparison to the mean-field value J? MF . Goldstone's theorem even requires the 
existence of a soft mode for the isotropic model we study here; its absence is the 
most serious failure of mean-field theory. Because of its neglect of collective mag- 
netization fluctuations, mean-field theory also always overestimates the critical 
temperature J2^|58|. 

2.3 New theory 

In this article we present a theory p8|]29| , |3"c| of DMS ferromagnetism that ac- 
counts for both finite itinerant-carrier spin splitting and dynamical correlations. 
It goes, thus, beyond the RKKY picture and mean-field theory and allows us 
to analyze the ordered state and its fundamental properties. The starting point 
for this new theory is an analysis of the system's collective excitations that we 
discuss in the following sections. 

3 Derivation of independent spin-wave theory 

In this section we use a path-integral formulation and derive, after integrating 
out the itinerant carriers, an effective action for the Mn spin system. We expand 
the action up to quadratic order to obtain the propagator for independent spin 
excitations. The latter will be the starting point for Section |J where we identify 
the system's elementary spin excitations and their dispersion. 
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3.1 Effective action 

Wc represent the impurity-spin density in terms of Holstein-Primakoff (HP) 
bosons p9| , 

S + (r) = ^2N M nS-tf(r)b(r)^j b(r) (3) 

S-(r) = 6t( P )^2JVMn5-fet(r)6(r) (4) 
S z (r) = N Mn S - b\r)b(r) (5) 

with bosonic fields b^(r), b(r). As can be seen from Eq. (||), the state with fully 
polarized Mn spins (along the z-direction, which we chose as the quantization 
axis) corresponds, in the HP boson language, to the vacuum with no bosons. 
The creation of an HP boson reduces the magnetic quantum number by one. The 
square roots in Eqs. ^ and ([|) restrict the boson Hilbert space to the physical 
subspace with at most 2S bosons per Mn spin. This constraint obviously yields 
an interaction between spin excitations. Later, when we go to the independent 
spin- wave theory, we will approximate the square roots by y/2NM n S, i.e., we will 
treat the spin excitations as free bosons. This is a good approximation as long 
as the spin-excitation density is small. 

We can write down the partition function of our model as a coherent-state 
path- integral in imaginary times, 

Z= (v[zz]V[W]e- So ( 6 ) 



with the Lagrangian L = Jd 3 r [zd T z + J2 a &<rd T & a ] + H(zz, WW). The bosonic 
(impurity spin) and fcrmionic (itinerant carrier) degrees of freedom are repre- 
sented by the complex and Grassmann number fields, z, z and W, W , respectively. 

Since the Hamiltonian is bilinear in fcrmionic fields, we can integrate out the 
itinerant carriers and arrive at an effective description in terms of the localized 
spin density only, Z = J T>[zz] exp(— [zz] ) with the effective action 

Ses[zz] = J dr Jd 3 r {zd T z — g/iBB(NMnS ~ zz)] 

-lndet [(G^)- 1 +8G- 1 {zz)\ . (7) 

Here, we have already split the total kernel G _1 into a mean-field part (G MF )~~ 1 , 
which does not depend on the fields z and z, and a fluctuating part 5G _1 , 



5G- 1 = ^ 
2 



(zt~ + zT+)^2N Mn S - zz- zzt 7 - , (9) 



where A = NMnJpdS — g* (J-b B is the zero-temperature spin-splitting gap for the 
itinerant carriers. We have allowed for the possibility of an external magnetic 
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field B along the ordered moment direction. Although the itinerant carriers have 
been integrated out, their physics is still embedded in the effective action of the 
magnetic ions. This fact is responsible for the retarded and non-local character 
of the interactions between magnetic ions. 

3.2 Independent spin-wave theory 

The independent spin-wave theory is obtained by expanding Eq. (Q) up to 
quadratic order in z and z, i.e., spin excitations are treated as noninteract- 
ing HP bosons. This is a good approximation at low temperatures, where the 
number of spin excitations per Mn site is small. 

We expand the term lndet [(G MF ) _1 + (5G _1 ] up to second order in 5G~ X , 

lndet (G- 1 ) = trln (G MF ) _1 +tr (G^G" 1 ) - ±tr (G MF SG^G^ SO' 1 ) , and 
collect all contributions up to quadratic order in z and z. The diagrammatic 
representation of the contribution linear and quadratic in <5G -1 are depicted 
in Fig. pi. We obtain (in the imaginary time Matsubara and coordinate Fourier 




Fig. 2. Diagrammatic representation of tr (G MF <5G _1 ) and tr (G mf 5G" 1 G mf <5G _1 ) . 
The solid lines represent the mean-field Green's function for the itinerant carriers and 
the incoming and outgoing dashed lines stand for z and z, respectively. 



representation) an action that is the sum of the temperature-dependent mean- 
field contribution and a fluctuation action, 

S e s[zz] = — z{k,V m )D~ l (k,l>m)z{k,V m )- (10) 

\h\ <fco ,m 

A Debye cutoff ko with fc^, = 6tt 2 Nmh ensures that we include the correct 
number of magnetic-ion degrees of freedom, |fe| < ko- The kernel of the quadratic 
action defines the inverse of the spin-wave propagator, 



D 1 (k,v m ) = -iv m + gfi B B 



-^gf^ E G f F (9> "n)Gf F (q + k, co n + u m ) , (11) 
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where G^f F (q, u) n ) is the mean-field itinerant carrier Green's function given by 
Gf F (q,uj n ) = -[iw n - (e q + aA/2- and e q = h 2 q 2 /{2m*). The frac- 

tional free-carrier spin polarization is £ = (p± — p^)/p. 



4 Elementary spin excitations 

We obtain the spectral density of the spin-fluctuation propagator by analytical 
continuation, iv m — > ft + i0 + and A(p, J?) = Im_D(p, J?)/tt. In the following we 
consider the case of zero external magnetic field, B = 0, and zero temperature, 
T = 0. We find three different types of spin excitations. 



4.1 Goldstone-mode spin waves 

Our model has a gapless Goldstone-mode branch reflecting the spontaneous 
breaking of spin-rotational symmetry. At long- wavelengths, the free-carrier and 
magnetic-ion spin-density orientation variation is identical in these modes. The 
dispersion of this low-energy mode for four different valence-band carrier con- 
centrations p is shown in Fig. (solid lines). At large momenta, k — > oo, the 
spin- wave energy approaches the mean-field result (short-dashed lines in Fig. ||) 

^Eoo = ^ MF = ^- (12) 

Note that the itinerant-carrier and magnetic-ion mean-field spin splittings differ 
by the ratio of the spin densities 

x=pi/{2N Mn S); (13) 

x is always much smaller than 1 in (III,Mn)V ferromagnets. Expansion of the 
T = propagator for small momenta yields for the collective modes dispersion, 

where is the Fermi energy of the majority-spin band. The spectral weight of 
these modes at zero momentum is 1/(1 — x). In strong and weak-coupling limits, 
A 3> ep and A <§C ep, respectively, Eq. ( O ) simplifies to 

Q k ] - -^£fc + 0(fc 4 ) for Z\»£F, (15) 

1 — x 

At long wavelengths, the magnon dispersion fit = pk 2 /M in an isotropic 
ferromagnet is proportional the spin stiffness p divided by the magnetization M, 
in our case M = NmhS — p£,/2. Note that the spin stiffness depends differently 
on the exchange constant J p d, effective mass m*, and the concentrations A^Mn 
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Fig. 3. Spin-wave dispersion (solid lines) for J p( j = 0.06eVnm 3 , m* = 0.5m e , Nmu = 
lnm" 3 , and four different itinerant-carrier concentrations p — 0.01 nm" 3 , 0.035 nm" 3 , 
0.1 nm" 3 , and 0.35 nm" 3 . The ratio A/e F is 2.79, 1.21, 0.67, and 0.35, which yields the 
fractional free-carrier spin polarization £ as 1, 1, 0.69, and 0.31. The short wavelength 
limit is the mean-field result fl M¥ = xA (short-dashed lines). For comparison, we show 
also the result obtained from an RKKY picture (long-dashed lines). 

and p in the strong and weak-coupling limits. The smaller the spin stiffness, 
the smaller the cost in energy to vary the spin orientation across the sample. A 
small spin stiffness will limit the temperature at which long-range order in the 
magnetization orientation, i.e. ferromagnetism, can be maintained. We see from 
Eq. ( |l5|) that the spin stiffness becomes small in systems with large effective 
band masses (such that A 3> £f), the same circumstances under which the 
mean-field-theory critical temperature becomes large. It follows that a mean-field 
description becomes more and more inappropriate |2^,[58| in systems with larger 
itinerant-carrier densities of states. This observation is important in devising 
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strategies for finding (III,Mn)V ferromagnets with larger critical temperatures. 
We will discuss this point in more detail in Section |^. 

4.2 Stoner continuum 

We find a continuum of Stoner spin-flip particle-hole excitations. They corre- 
spond to flipping a single spin in the itinerant-carrier system and, since a; < 1, 
occur at much larger energies near the itinerant-carrier spin-splitting gap A 
(see Fig. ||). For A > ep and zero temperature, all these excitations carry spin 
S z = +1, i.e., increase the spin polarization. They therefore turn up at negative 
frequencies in the boson propagator we study. (When A < ep, excitations with 
both S z = +1 and S z = — 1 contribute to the spectral function.) This contin- 
uum lies between the curves — A — e k ± ^\f^k^F and for A < tp also between 
-A + e k ±2^/e k {e F -A). 



1.05 




0.95 1 ' 1 1 1 1 1 ' 1 1 1 

0.000 0.001 0.002 0.003 0.004 0.005 

k/k D 



Fig. 4. Stoner excitations and optical spin-wave mode in the free-carrier system for 
J pd = 0.06eVnm 3 , m* = 0.5m e , iV M n = lnm -3 , and p = 0.35nm -3 . In an RKKY 
picture these modes are absent. 



4.3 Optical spin waves 

We find additional collective modes analogous to the optical spin waves in a 
ferrimagnet. Their dispersion lies below the Stoner continuum (see Fig. ||). At 
small momenta the dispersion is 

-w-w-i-^fa-'-ds^wyw. (.7, 
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The spectral weight of these modes is —x/(l — x) at zero momentum. 

The finite spectral weight at negative energies indicates that, because of 
quantum fluctuations, the ground state is not fully spin polarized. 

5 Comparison to RKKY picture 

For comparison we evaluate the T = magnon dispersion assuming an RKKY 
interaction between magnetic ions. This approximation results from our theory 
if we neglect the spin polarization in the itinerant carriers and evaluate the 
static limit of the resulting spin- wave propagator defined in Eq. (|ll]). The Stoner 
excitations and optical spin waves shown in Fig. ^ are then not present and the 
Goldstone-mode dispersion is incorrect except when A <C ep, as depicted in 
Fig. H (long-dashed lines). As a conclusion, the RKKY picture can be applied 
in the weak-coupling regime only (in this limit the long-wavelength fluctuations 
are described by Eq. (|l6|)). Outside that regime it completely fails as a theory 
of the ferromagnetic state. 

6 Comparison to mean-field picture 

In the mean-field picture, correlations among the Mn spins are neglected. The 
mean- field theory can be obtained in our approach by taking the Ising limit, 
i.e., replacing S ■ s by S z s z . This amounts to dropping the last term in Eq. ( [Tl|) 
or the second diagram in Fig. [| It is this term that describes the response of 
the free-carrier system to changes in the magnetic-ion configuration. When it is 
neglected, the energy of an impurity-spin excitation is dispersionless, Q Ml = xA, 
and always larger than the real spin-wave energy as can be seen in Fig. ^ (the 
short dashed line is the mean-field value and the solid line shows the real spin- 
wave energy). 

However, if both the weak-coupling (A/ep <C 1) and dilute-density (p/NMn <C 
1) limit are met, the spin- wave energy is almost dispersionless Slf. j? MF (see 
Fig. ^|), except in the narrow window near k — which is protected by Gold- 
stone's theorem. In this regime, the spatial correlations among Mn spins are less 
important. 

7 Comparison to a ferrimagnet 

Some features of the excitation spectrum are, not coincidentally, like those of a 
localized spin ferrimagnet with antiferromagnetically coupled large spin S (cor- 
responding to the magnetic ions) and small spin s (corresponding to the itinerant 
carriers) subsystems on a bipartite cubic lattice (see Fig. |^). We can represent 
the spins by HP bosons and expand up to quadratic order. Then we either di- 
agonalize the resulting Hamiltonian directly using a Bogoliubov transformation 
or, as in our ferromagnetic-semiconductor calculations, integrate out the smaller 
spins using a path-integral formulation. The latter approach is less natural and 
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Fig. 5. Spin-wave dispersion (solid line) for A/cf = 1/3 and p/Nmu = 0.001, in 
comparison to the mean-field value O mf (short-dashed line) . 

unnecessary in the simple lattice model case, but leads to equivalent results. In 
both cases we find that there are two collective modes with dispersion 



In analogy to the ferromagnetic DMS, we defined the ratio of the spins as 
x = s/S, the mean- field spin-splitting gap is A = 6JS for the smaller spins, 
where J is the exchange coupling, and 7*. = (1/3) JZJl — cos(kia)] with lattice 
constant a. The two collective modes correspond to coupled spin waves of the 
two subsystems. One is gapless, the other one gapped with A(l — x), exactly 
as for ferromagnetic semiconductors. The bandwidth is xA, and the spectral 
weights at zero momentum are 1/(1 — x) and —x/(l — x), respectively, as in our 
model. 



Fig. 6. Localized spin ferrimagnet. Local moments with large spin S and small spin s 
are sitting on a bipartite cubic lattice. They are coupled by a nearest-neighbor antifer- 
romagnetic interaction. Only one dimension is shown. 



4 1)/(2) - y H 1 _ x) ± y/(i-x) 2 + te-y. 



(18) 
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This comparison demonstrates that some of our results for the collective 
excitations of ferromagnetic semiconductors do not depend on details of the 
model but reflect generic features of a system in which two different types of 
spins are antiferromagnetically coupled. 

8 Alternative derivation of the spin-wave dispersion 

The energy of long-wavelength spin waves is determined by a competition be- 
tween exchange and kinetic energies. To understand this in detail we impose the 
spin configuration of a static spin wave with wavevector k on the Mn spin sys- 
tem (see Fig. [?]), evaluate the ground-state energy of the itinerant-carrier system 
in the presence of this exchange field, and compare this with the ground-state 
energy of a uniformly polarized system. We will rederive the dispersion given 
in Eq. ( |l4| ) as the ratio of energy increase SE and decrease of the total spin SS. 



\ \ 

v — 1 — i 




Fig. 7. a) To minimize the kinetic energy the itinerant carriers align antiparallel to the 
average Mn spin orientation. In this case, exchange energy is increased, b) To minimize 
the exchange energy the itinerant carriers align antiparallel to the Mn spin orientation 
everywhere in space. In this case, kinetic energy is increased. 
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How will the itinerant carriers respond to the inhomogeneous Mn spin ori- 
entation? On the one hand, they can minimize their kinetic energy by forming 
a homogeneously polarized state with the spin pointing in negative z-direction 
(see Fig. 0a). In this case, however, the Mn and the itinerant-carrier spins are 
not aligned antiparallcl, i.e., the exchange energy is not optimized. One the other 
hand, if the valence-band spins decide to follow the Mn spin orientation every- 
where in space, as indicated in Fig. ^Jd, then the exchange energy is minimized 
but the kinetic energy is increased because of a new contribution to spatial vari- 
ation in the itinerant-carrier spinors. Since the energy scales of the exchange and 
kinetic energy are provided by A and ep, respectively, it is the ratio A/ep that 
determines which term controls the itinerant-carrier behavior. 



8.1 Strong-coupling limit 

In the strong-coupling regime, A/ep 3> 1, the itinerant-carrier spins minimize 
the exchange energy by following the Mn spin orientation in space as shown in 
Fig. Qd. The energy SE will therefore be purely kinetic and the spin- wave energy 
will not depend on A in this regime. To evaluate the kinetic energy cost we use 
a variational ansatz with trial single-particle spinors 

for a free carrier with wavevector q, where k is the wavevector of the spin wave 
going in ^-direction, and 9 is the angle by which the Mn spins are tilted from 
their mean- field orientation. The total kinetic energy, i.e., the sum of the single- 
particle kinetic energies £'kin,a('Z) = (^(<z)| — ?i 2 V 2 /(2m*)|!f r (q)) a , is minimized 
for a m j n = (k/2) cos 6, which yields for the increase in kinetic energy SE due to 
the spin wave 

SE 1 

— = -pe k sin 2 8 . (20) 

We note that 

y = (N Mn S-p/2)(l-c 0fi 8) (21) 

and obtain 

fl« = ^ £feC os 2 (fl/2). (22) 
1 — x 

In the limit of small spin-excitation density, 9 — > 0, we recover the independent 
spin- wave theory result, Eq. (|l5|). 



8.2 General case 



We now generalize this discussion to arbitrary values of A/ep. The Hamiltonian 
for itinerant carriers in an effective potential generated by the Mn spins is given 



Theory of Ferromagnetism in Diluted Magnetic Semiconductors 15 



in second quantization by 
q 

«,± 

with quasiparticle energies 



.± + (4/2) cost 
2 (4/2)sin(9 



(4/2) sin6> 
e q+# -(4/2)cos( 



(23) 



q=p tt cos £ 



77 sin C7 



4 sin 2 6> /?i 2 (fc • q) 
~2 + 4 V 2^ 



0(fc 4 ). (24) 



It is now straightforward to calculate the band energy of the new quasiparticle 
bands. In the strong-coupling limit, only the lower quasiparticle band is occupied, 
the only term from Eq. (pi]) which enters the energy difference SE is et sin g, and 
we recover Eq. (pcj) . For arbitrary values of A/ep we arrive after a lengthy 
calculation at 



SE 



1 



— = —pet sin 



which yields together with 



y = (7V M nS-K/2)(l-cos< 



the desired result (compare to Eq. (14)) 



i'k — i £ fc 

1 — a; 



(25) 



(26) 



(27) 



9 Extensions of the model 



The discussion in preceding sections has been based on a number of assumptions 
and idealizations that have enabled us to expose some essential aspects of the 
physics in a relatively simple way. Among these we mention first the assump- 
tion that Mn ion charge fluctuations are weak and therefore can be described 
perturbatively, representing their influence by effective exchange interactions 
with itinerant carriers. Failure of this assumption would completely invalidate 
the approach taken here. In our judgement, however, experimental evidence is 
completely unambiguous on this issue and favors the approach we have taken. 
Another essential assumption we have made is that we can use an effective-mass 
approximation, or more generally an envelope-function approach, to describe 
the influence of the kinetic-exchange interaction on the itinerant bands. This as- 
sumption is safe provided that itinerant-carrier wavefunctions are not distorted 
on atomic length scales by their interactions with the Mn acceptors. We believe 
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that effective-mass-theory's success |Q in interpreting isolated Mn acceptor 
binding energy data, albeit with central cell corrections (see below) not included 
in the simplest model, establishes its applicability in (III,Mn)V semiconductors. 
Many other extensions of our simplest model that may be required for experi- 
mental realism can be included without any essential change in the structure of 
the theory. We mention some of these below. 

• The valence bands are p-type, with 6 orbitals including the spin-degree of 
freedom. In addition, they are not isotropic, but reflect the cubic symmetry 
of the crystal. A more realistic description is provided by a six-band Kohn- 



Luttinger Hamiltonian 1 40 that includes an essential spin-orbit coupling 
term and is parametrized by phenomenological constants whose values have 
been accurately determined by experiment. Furthermore, depending on the 
substrate on which the sample is grown, either compressive or tensile strain 
modifies the band structure. This effect can also be included in the Kohn- 
Luttinger Hamiltonian. 

Mean-field calculations based on this more realistic band structure have been 



performed recently 34 36J ■ Since the model is no longer isotropic the spin- 
wave dispersion will be gapped. 

The impurity density is not continuous. Instead, discrete Mn ions are ran- 
domly distributed on cation lattice sites. In recent numerical finite-size cal- 
culations we have found |38| that this disorder has important effects. It leads, 
e.g., to an increase of the spin stiffness in the strong-coupling regime and 
presumably also plays an important role in determining the T = resistivity. 
Evidence from studies of (II,Mn)VI semiconductors suggests that direct an- 
tiferromagnetic exchange between Mn ions is important when they are on 
neighboring sites. The way in which this interaction arises from the micro- 
scopic atomic length scale physics is incompletely understood. This interac- 
tion will become more and more important for increasing doping concentra- 
tion. 

Coulomb interactions between free carriers and Mn acceptors and among the 
Mn acceptors are also important. For dilute Mn systems, the free carriers 
will reside in acceptor-bound impurity bands. In ignoring these interactions 
in our simplest model we are appealing to the metallic nature of the free 
carrier system at larger Mn concentrations, that will result in screening. We 
believe that this source of disorder will be important in determining the 
T = resistivity and the magnetic anisotropy energy. 

Including Coulomb interaction among the itinerant carrier enhances ferro- 



magnetism as shown in mean-field calculations 1 33 . 

Microscopic details of the p-d exchange physics suggest a finite effective 
range of the exchange interaction pif . In fact, a finite range is also required 
as an ultraviolet cutoff since a delta- function type of interaction with discrete 
Mn spins leads to an ill-defined problem with infinite negative ground-state 
energy. In this article we have circumvent this problem by assuming a ho- 
mogeneous Mn density. 
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To estimate the additional effect of a finite exchange interaction range we 
replace the delta-function with a Gaussian function 

J -<'> ^ (2^W^ 2/(2 ' 2) < 28 > 

with cutoff parameter I. We can redo the analysis of the spin waves and 
find that, in the long- wavelength limit, the dispersion is modified from Qk = 
Dk 2 +C(fc 4 ) to 

n k = (D + f2 MF l 2 )k 2 +0(k 4 ), (29) 
i.e., the spin stiffness acquires an additive constant. 

10 Conclusion 

In conclusion we have presented some aspects of a theory a ferromagnetism in 
doped diluted magnetic semiconductors which goes beyond both the RKKY pic- 
ture and beyond the mean-field approximation for the kinetic-exchange model. 
The present discussion has focussed on the dependence of the system's collec- 
tive spin excitations on model parameters and on implications for long-range 
magnetic order. 
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